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We consider a quasi-linear Hamiltonian system with one and a half degrees of 
freedom. The Hamiltonian of this system differs by a small, ~ e, perturbing term 
from the Hamiltonian of a linear oscillatory system. We consider passage through a 
resonance: the frequency of the latter system slowly changes with time and passes 
through 0. The speed of this passage is of order of e. We provide asymptotic formulas 
that describe effects of passage through a resonance with an accuracy 0(ea). This is an 
improvement of known results by Chirikov (1959), Kevorkian (1971, 1974) and Bosley 
(1996). The problem under consideration is a model problem that describes passage 
through an isolated resonance in multi- frequency quasi- linear Hamiltonian systems. 



1 Introduction 

Study of passage through an isolated resonance in a multi-frequency quasi-linear Hamiltonian 
system can be reduced to the case of one- frequency system (see, e.g., [1]). The corresponding 
Hamiltonian has the form 
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Here J, tp mod 2tt are conjugate canonical variables, r is a slow time, f = e, and e is a small 
parameter, 0<£<1. Equations of motion are 



j = _ c BH ! ^f) )+e 9g,(/,V,r) 

For e = 0we get an unperturbed system with the Hamiltonian Hq(I, t) = oj(t)I and action- 
angle variables /, tp. The function u is the frequency of the unperturbed motion. For some 
value of the slow time r*, where there is a resonance, u vanishes: cj(t*) = 0. We assume that 
the resonance is non-degenerate: u'^ = w'(r*) 7^ 0. Here "prime" denotes the derivative with 
respect to r. Let, for definiteness, u[ > 0. We assume that r* is the only resonant moment 
of the slow time: oj(t) is different from Oatr^r,. 

Action I is an adiabatic invariant: its changes along trajectory of (1.2) are small over long 
time intervals. For motion far from the resonance value I oscillates with an amplitude ~ e. 
Passage through a narrow neighbourhood of the resonance leads to a change in I of order 
\fe (so called jump of the adiabatic invariant). There is an asymptotic formula for this jump 
([1], [2]). Let I_ and J + be values of I along a trajectory of (1.2) at moments of slow time 
t_ and t+, where r_ < t* < r+. Then 

+00 / 

7+-/_ = -Ve J d6 + 0(e). (1.3) 

—00 

Here is the value of tp on the considered trajectory at r = r*. There are formulas for 
change of the angle (phase) tp due to passage through the resonance as well [1]. 

One can replace in the left hand side of (1.3) I± with values of the improved adiabatic 
invariant J±, but the error estimate in (1.3) still will be 0(e). It was suggested in [3] to 
eliminate an asymmetry in (1.3) by replacing in the right hand side J_ with where is 
the value of I on the considered trajectory at r = r*. A numerical simulation in [3] shows 
that this symmetrization indeed improves the accuracy of formula (1.3) for I± replaced with 
J± considerably. It is conjectured in [3] on the basis of the numerical simulation that the 
error term in the modified formula is 0(e^). 

In the current paper we prove this conjecture by means of a Hamiltonian adiabatic per- 
turbation theory. We show that this improvement of accuracy occurs due to cancellations 
of many terms in formulas of the perturbation theory considered up to terms of 4th order 
in e. We obtain also formulas which describe change of phase tp due to passage through 
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resonance with the same accuracy 0(e^). As a result, we obtain formulas which allow to 

3 

predict motion in post-resonance region with accuracy 0(e^), provided that the motion in 
the pre-resonance region is known. In the last section we provide a numerical verification of 
these formulas. 
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2 Main theorems 



We consider Hamiltonian system (1.2) with Hamilton's function (1.1). We assume that the 
function H is of class C 4 for (J, ip, r) 6 D = Dj x R x D T , where Dj and D T are some open 
intervals in R. We assume that Hi is 27r-periodic in ip and that the frequency u does not 
vanish in D T other than at r = r*. At the resonance state, cu(r*) = 0, but u/(r*) = oj'^ ^ 0. 

Let I(i), y>(£) be a solution of (1.2) on a time interval [£_,i + ], where t± = t±/e, and t± are 
some constants. Let r_ < r* < r + . We denote J± = I(t±), <p± = <p(t±), and i* = I(t*), 
<p* = ip(t*). Here t* = t*/e. 

Theorem 2.1. 

+ 00 / 

J + + £Mi(/ + ,^ + ,r + ) = J_ +£%(/_, v?-,r_) - y d^ + 0(£2), 



(2.4) 



-oo 



n 1 7 / ^ ?9i7 1 (J_,r) J T fdH 1 (J+,T) J 
p + + ev 1 {I + ,cp + ,r + ) = <p- + ev!(I-,<p-,T-) + - J u{r)dT+ / 0/^ 



T— 



^ ' ar d^ + p.v. £ y — — — dr 

ellneS 3 ^ 2 )^ 1 *^*) T^i( J *> ^ + ^g!l<> d g , Q(g |) (2 5 ) 



4c< <9/ 3 J <9^ 



Theorem 2.2. 



1 _ f dH x {h,^ + ^d 2 



T* 



I± ± 2 ""^"^ 2 ' d ^ + (2-6) 



-oo 

r* r* _ +oo 



V?* = (p± + - I w(r) dr + / — dr =f ->/e / d ^ 
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+ ^«^) +0(£) (2?) 
In the above theorems, 

u x = — , Jfi = i?i - ffi, #i = (.Hi) , vi = - 7^ dip+ ( — dV> 
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2tt 

and angular brackets denote averaging with respect to ip: (f)^ — w- J /(v 9 ) d</?- 

71 o 

Moreover, 



J± = I± + eu 1± , n 2 {J±,T) = - — — — , U 1± = U 1 {I±,Lp±,T±). 

2a; (r) dl z 

The new results here are: 

• estimate 0(e%) in the first formula of Theorem 2.1, 

• the last two terms and the estimate 0(e*) in the second formula in Theorem 2.1, 

• the last term and the estimate 0(e) in the second formula in Theorem 2.2. 

Combining results of Theorems 2.1, 2.2, we can obtain prediction of motion in the post- 
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resonance region with accuracy 0(e%) as follows: 
Corollary 2.1. 

+oo / 

/+ = /_+eu 1 (/_,<^_,r_)-eui(/_,^ + ,r + )- ^fej ^ — ^— ^ J -d8 + 0(e*), 

—oo 

tr ~ \ l l ,m 1 : 9Hi(J-,r) j 7dH 1 {J + , 
tp+ = (p-+ev 1 (I-,tp-,T-)-ev 1 (I-,(p + ,T+) + - o;(r)dr+ / — dr + / ^— 

+ 0O ~ „ / 

dH^^ + ^-e 2 ^) Aa [dK 2 (I-,r 



T— 7* 



d6> + p.v. e / — dr 



dl J dl 

T— 

el hxe&ffiYil-, n) TdHtil^, ft + ^# 2 , r, 



+oo 

d# + 0(^), 
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where 



0+ = ^_ + - / u;(r)dr + / dr, 



+00 / 
1 /■ <9#i(/_,</3* + ^ 2 ,r*) 



/* = I„--Ve / ^ r ; dfl, 

2 y ^ 



-00 



t» r, _ +00 ~ „ , 

= <P- + -J u(t) dr + J — dr + - J ^—t d6 

T— T— — OO 

g ln g (9 2 (gy(I^) 
+ Aui dP 

1 ? / n j ?0#i(/-,t) , 
y>* = </?- + "/ ^( r ) dr + / — dr, 

+OO v , 



3 Procedure of adiabatic perturbation theory 

In order to get the above estimates, four steps of adiabatic perturbation theory are performed 
in both original Hamiltonian system (in subsection 3.1) and an approximate Hamiltonian 
system (in subsection 3.2). We follow an approach of [4] here. 

3.1 Original Hamiltonian system 

Consider dynamics described by the Hamiltonian 

H{I,<p,T)=L>(T)I + eH 1 {I,< P ,T). 

The frequency oj{t) can be expanded near the resonance as 
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w(r) = ui(r - r.) + ~w?(r - r*) 2 + 0(r - r*) 3 



with u'^ ^ 0. 
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So we can get formula for if as 

T 

V^ = ^ + ^(^:(r-r,) 2 + ic:'(r-r,) 3 + 0(r-r^+ J dH ^ T ^ dn . ( 3 .10) 
Let 

eS{J, if, r, e) = eS x {J, if, r) + e 2 S 2 {J, if, r) + e 3 S 3 (J, if, r) + e*S 4 (J, if, r), (3.11) 
where S is 27r-periodic in f and (S)^ = 0. Here we have not defined Sj yet. 

Make the canonical transformation of variables (J, 99) 1— >■ (J, ^) with the generating function 
J f + eS(J, f, t, e). Old and new variables are related as follows: 

I = J + e dS{J >?' T > £ \ ^ = f + e dS{J t T ^\ (3.12) 
dip oJ 

The new Hamiltonian, which describes dynamics of variables J, ip, is 

/ dS\ dS OS 

n(J,if(J,ip,r),r) = u(t) [J + £—) +eH!(J + £—,if,r) +s— 

. .dS TT . T , 2 dH x {J,if,T)dS 2 dS 
+ e^d^H 1 {J 1 ip 1 T) (9S\ 2 i e ^H^J^r) fdS_V 



2 dP \dipj 3! a/ 3 \dipj 

We would like to find S^, j = 1, ...,4 such that there is no dependence on the new phase if) 
in the new Hamiltonian T-L in terms up to 4th order in e. Thus the new Hamiltonian should 
have a form 

H{J, f (J, if), t),t)= u{t)J + e1l 1 (J, r) + e 2 TZ 2 { J, r) + e 3 Tl 3 (J, r) + e 4 Jl i (J, r) + e 5 TZ 5 (J, ip, r). 

(3.14) 

Here we have not defined functions 71 j yet. 

OS- 

Equating terms of the same order in e in (3.13) and (3.14) we find that functions — — have 

Of 

forms: 
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dS\ a 1 (J, <p,r) 



dip 



dip 



dip 



dip w(t) 

dS 2 uj\r)a 2 {J^,r) , (3 2 (J,ip,r) 



u 3 {r) 

,2 



dS 3 _ {c/(T)yaP(J,<p,Tl + p*(J,<p,T) 



w 5 (r) 

,3 



dS, _ ( W '(r))V(J,^r) + /3 4 (J,^r) 



uA(r) 



where a , cr, a , a , f3 , /3 , /3 are smooth functions. 
Denote 

f(J,<p,r)=J f(J,<p,r) d<p-(^J f(J,<p,r) dip 
We can find the explicit form of 5^: 

1 



51 = — Hi, 

UJ 

52 = nHl + 



^4 = 7~ -"1 



Co' 



Co' 
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satisfying (S^ = 0, j = 1, 2, 3, 4. 
So we have the expressions 

dSi e dH 1 
dJ u dJ 



-2 ap2 



S 



dJ 
dJ 



4^4 
: dJ 



co 3 dJ + 


u 2 dJ ' 


3uj' 2 dHi 


e 3 dp 3 


UJ 5 dJ 


f <9J ' 


15a;' 3 dHi 


e 4 <9/3 4 



cj 7 <9J dJ ' J 
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Also 



K 1 {J,t) = (H 1 {J,< P ,t)) 

dH 1 (J,<p,T) dS 1 V p 
dl dtp I 

ld 2 H 1 {J^,r) fdSA 2 dH 1 (J,ip,r)dS t X 



TZ 2 (J,t 
TZ 3 {J,t 

tz 4 (j,t 



2 dP V dtp J dl dip 

/l ggxQ^r) /dSA 3 | d 2 H 1 (J,<p,r) dS 1 dS 2 | dH 1 (J,ip,r) dS 3 X * 
\ 6 dP \dip J dP dip dip dl dip 



The new Hamiltonian is 

H(J, <p(J, ip, t),t) = uj(t)J + eKx + e 2 K 2 + e 3 K 3 + e 4 TZ 4 



e 5 0(l) , 5 dS 4 



uj 7 (t) dr 

For |r — r*| > v^e, the motion is described by differential equations: 
dU dip &H fdtpy 1 



J 



dip dip dip \dip 



e 



dip \ w(r) w 3 (r) w 5 (r) w7 ( r ) 

B 7(a/(r))V(J,^,r) 57 (J,^,r)\ / A kl k (J,tp,r) 



-i 



C0 8 (t) L0 7 (t) j\ 0J 2k ~ l (T) \{t-T^ 



s 2 



■ _m mdip_m &h (_ (fs\ ( d 2 s y 1 

^~'dJ + 'd^'dJ~'dJ + 'dlp~' \ £ dPJ\ +£ dmp~) 

dH dH d ( dS 1 +£ 2dS 1 + £3 dS l + £ ,dSA( 1 + fe 



dJ dip dJ V 9J dJ dJ dJ J \ \tu 
9 ( „ 2 ^ 4^ e 5 0(l) .dS 4 



V ^V) ^ 

^ + *X\ (X + + (l + ol 



UJ a LO 7 J \ LO UJ 3 OO 5 LO 7 J \ \0J 

u(r) + A fa + £ *n 2 + e*K 3 + e*K 4 + e^) + £ + O ( 



(3.18) 



.(t-t*) 1( V 
(3.19) 

Here 7, 7 1 , . . . , 7 9 , 7 1 , . . . , 7 4 are smooth functions. 



S 



Remark. By differentiating if) = tp + e ' ^ - — — — — - with respect to J on both sides, 

oJ 

dip d 2 S d 2 S dtp m „ d 2 S ( d 2 S \ x 

we obtain = - + e— 2 + Sgj^gj- Therefore, - = -e^ ■ 



3.2 Approximate Hamiltonian system 

Now let us consider the approximate Hamiltonian 

H(I a , (p a , t) = Q(r)I a + eHi(h, ip a , r*) 

Equations of motion are 

■ dHi(h,<p a ,r*) 
4 = -e , <p a = U{t). 

Here Q(r) = u'^r - r*), ^ ^ 0. 

We will consider the solution of these equations with initial conditions at resonance, i.e. 
when r = r*: I a (n) = I* = I(U), ip a (T*) = y?* = <p(U). We get the formula for cp a as 

y? a = + ^ UJ '*( r - T *f- ( 3 - 20 ) 

Let 

eS a (^ r, 5) = eSZ(tp a , r) + £ 2 S 2 a (<^, r) + £ 3 S 3 > a , r) + e 4 ^(^ a , r) (3.21) 
where S" 1 is 27r-periodic in (p a and (S a ) Va = 0. Here we have not defined yet. 

Make the canonical transformation of variables (J , <y2 a ) i— >■ ( J a , ip a ) with a generating function 
J a fa + £S a ((p a , t, e). The old and new variables are related as follows: 

4 = Ja + ^ ~ , = fa- (3.22) 

The new Hamiltonian, which describes dynamics of variables J a , (p a , is 
n a (J a ,ip a ,r) = £l(<r)[j a + e-—)+£Hi(h,<p a ,n) + £- 



dS a dS a 
= fl(T)J a + eQ(T)— + eH 1 (h,ip a ,T*)+e 2 —. (3.23) 
dip a Or 

We would like to find Sj, j = 1, ...,4 such that there is no dependence on the phase y? a in 
the new Hamiltonian T-L a in terms up to 4th order in e. Thus the new Hamiltonian should 
have a form 

H a (J a , <p a , t) = fi(r) Ja+eKKJa, r)+e 2 TZ a 2 ( J a , r)+e 3 TZ a 3 (J a , r)+e i Ul(J a , T)+e 5 K a 5 (J a , <p a , r). 

(3.24) 
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dS a 

Equating terms of the same order in e in (3.23) and (3.24) we find that functions — — have 
forms 





1 / T \ 

a L {I*,(p a ,n) 




fi(r) ' 


ds a 2 


Q'(T)a 2 (I*,(p a ,n) 




fi 3 (r) 


OS" 


(Q'(r)) 2 a 3 (/*,^ a ,r*) 






dSl 









The new Hamiltonian is 

U a (J a , <p a , t) = Q(r)J a + e{H x {h, <p a , n)f a + e 
The motion is described by differential equations 



<ta , 5 dS2{<pa,r) 



dr 



dip a dT 

k a = fi(r). 
Here 7 is a smooth function. 



5 &SS(<p a ,T) _ 5 7(1T(t)) a*(I*,<p a ,n) | 5 7(/„^,r,) 



fi 8 (r) 



(3.25) 



(3.26) 



4 Proofs of the theorems 



We will prove asymptotic formulas for the action variable I in both Theorems 2.1 and 2.2 first, 
and then asymptotic formulas for the angle variable (p in these Theorems. Denote Ti = t* — 
5 1 / 2 , r r = t* ti jT = Ti^/e 1 / 2 . Denote u = — |^ and u a = — f^- In the following text, we 

use the notations Zj = Z(tj), where Z = I, I a , J, J a , <p, <p a , ip, S, S a , Sk, S%, u, u a , Uk, u%, v, Vk 
and tj = t-,ti, t*, t r , t + , k — 1, 2, 3, 4. 

For simplicity of the exposition we will assume that r + and r_ are symmetric with respect 
to t*: r + — t* = r* — r_ . General case can be easily reduced to this one. 



4.1 Principal lemmas 

The following lemmas will be used in the proof. 
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Lemma 4.1. 

(tu'(T)) a _ (n'(r)) a 0(1) const 0(1) 



u b (r) Q b (r) (r-r,)^ 1 (r - r*) b (r-r,) 6 " 1 ' 

where a,b G N. 
Lemma 4.2. 

I z -/* = 0(Vi), I r -I. = 0(^), I,-I r = 0(VS). 
J, -J, = 0(Vi), J r -h = 0{y/e), J,-J r = 0(Vi). 

Lemma 4.3. 

(Jj - J») + (J r - J,) = 0(e), (tpi - </?aJ + (<£V - far) = °{ £ )- 
Proof. • For [(J/ — J*) + (J r — /*)], substitute with formulas (3.12) and (3.15): 



dtp J \ dtp 

t, 

dH x (I, tp, t) ^ _ ^ a l (J u tp U Ti) _ ^ 2 a 2 {J h tp h Ti) _ ^ 3 a 3 (J h tp h Ti) _ ^ a 4 (J h tp h ri) 



dtp n - (ti - r*) 3 (n - r*) 5 (n - r*) 7 

dH x {I,tp,r) a?-(J T ,tp r ,T r ) 2 a 2 (J r ,tp r ,T r ) 3 a 3 (J r ,tp r ,T r ) A a 4 (J r ,tp rj T r ) 



-£ / dt-e £ — ; S — -. E — n= hO £ 

dtp r r - r* (r r - rj 3 (r r - r*) 5 (r r - r*) 7 



T* 



dH x (I^p a ,n) dr _ j dH x (I^p a ,r*) dr _ c k & k (Jh<Pi,Ti) _ ^ Q k & k (Jh<Pi,i~i) q, n 

= 0(5). 

For Utpi — tp ai ) + [tp r — tp ar )\, substitute with formulas (3.10) and (3.20): 

(tpi - tp ai ) + (Pr ~ Pa r ) 

\* i / n31 , l ni f 9H x (I,tp,r x ) f dH x (I,tp,r x ) 

= ^ ifa - r *) + ( r r - r *) j + -ofa 7 *) + y — ^ — dr i + y — of — dr i 

7* T* 

dH x (I*,tp a ,n) [ dH x (I*,p a ,T*) 

dr x + / — dr x + O(e) 



9/ y dl 

O(e) 

□ 
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Lemma 4.4. 



e aK -—^f (I,p,r) = O(e a+ ^), 

U C [T j 



if t G {ti, T r }. Here a 6l, b,c G Z, and /o is a smooth function. 
The same estimate is valid if u(t) is replaced with Q(r). 



Lemma 4.5. 

t 



J co^eh) v 



if t E \t-,tj\. idere a G R, c G Z, c > 2, and /o is a smooth function. The same estimate is 
valid if ti is replaced with t r and t G [£ r ,£+]. 

The same estimate is valid ifu(r) is replaced with Q(t). 
Lemma 4.6 (Cancellation lemma near the resonance). 

(r - r*) 2 ^ 1 /!(/*, <p a , r*) dr = 
n 

where k G N and f\ is a smooth function. 

Lemma 4.7 (Cancellation lemma far from the resonance on symmetric intervals). 

fir - rn(I., „„r.) dr + /(r - r.)~ A </.,„., r.) d T = 
where k G Z and /2 is a smooth function. 

Lemma 4.8. Let f = f(I, tp, r) be a twice continuously differentiable function. Then 
f(I,<P,T) ~ f(I*,<Pa,T*) 

T / T 

df(I*,ip a ,n) f dHi(h,<p a ,T*) df(I*,<p a ,n) ( u"{t - r*) 3 f dH^I*, (p a , t* 

clri H h 1 



9/ J dcp dtp \ Ge J 81 

t* \ r* 

a/(J.,^,r.) (T _ t>) + 0(t _ r<)2 + (<I^)1) + o ((I^Z 
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Proof. 



df(I*,<Pa,T*), T r s , df(I*,ip a ,T*) df(I*, ip a , T») , . 

97 (/ " J * } + di (Lf ~ ^ + d~r (T " ^ 

+ 0(1 - hf + 0{ip - (fa) 2 + 0(r - r*) 2 

r / r 



9/ J dip dip \ 6e J dl 

+ ^(r-r t ) + 0(r-^ + 0f^V0 Ar "^ 



-2 



<9r \ e y \ £ 

^/(^*, </>a, T*) / OHi(h, ip a , T*) ^ + df(h,ip a ,n) I U)"(t - T*) 3 + /" 8Hi(I*, ip a , T*) ^ 



dl J dip dip \ 6e J dl 

+ 9/( V' ^ (r ~ r.) + Q(r - r,) 2 + O f +0^^ 



dr \ e J \ e 



-2 



□ 



Lemma 4.9. Let / = f(J,ip,r) be a twice continuously differentiable function. Then for 
r e [r-,Ti] U [r r ,T + ], 



f(J,<f,r) - f(h,ip a ,n) 

/ T 

df(h,ip a ,n) I f dH^I^ipa,^] 



dn 



81 \ J dip 

— d~p — — — + y — — dri + — ^ — (r - 



+ 0(e) + 0(r - r*) 2 + ^ ^- + O 



:r-rj 4 \ /(t-t*) 6 



5 2 



J X {U,Wa,T*) g2 f 2 {h,ip a ,n) c3 f 3 (I*,<Pa,n) ^ 4 f 4 (h,ip a ,T*) 

r - r* (r - r*) 3 (r - r*) 5 (r - r*) 7 

where f 1 , f 2 , f 3 , f A are smooth functions. 
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Proof. 



f(J,ip,r) - f{I*,ip a ,n) 

df{I*,<p a ,n) df(h,ip a ,r*) df(I*,<p a ,n) 

dl (J - Q + dj> ^ - ^ + d~T (T " 

+ 0(J - hf + 0(<p - ifaf + 0(t - r*) 2 

df(I*,<p a ,T*) ( dS(J,tp,r)\ df(h,ip a ,r*) df(I*,(p a ,n) 

i - 1* - e H [tp - (f a ) H (r - r* 



dl \ dip J dip dr 

+ 0(J - hf + 0{ip - ip a f + 0(t - r*) 2 

df(I*,<p a ,n) ( _ _ dS(J,tp,r) \ df{h,ip a ,T*) I u"(t - r,) 3 j dHijl^ip^n) 

dl { * £ dip J dip 1 6e J dl 

+ df(i,^,T, ) {T _ t>) + Q(e) + 0(t _ r>)2 + Q p^z£) + o ((iz_^ 6 ' 

Then 

r r 

J dip J dp ^ ~ T *' 



Also 

dS(J,ip,r) dS(I*,<p a ,T) 



0(e) + 0(r - r* 



<9<£> 

f\I*,<Pa,T*) 2 f 2 ( I *, l Pa,n) 3 f 3 (h,p a ,T*) 4 / 4 (/*,^a,n) 

= "^^ 6 (r-r,)* (r-r,)* (r - r.)? + ^ + 0(r " 

This implies the result of the lemma. □ 



4.2 Proof of formulas for action variable 
4.2.1 Principal identity for formula (2.4) 

Let 

t+ tr 

ei = J (J- j a ) dt, e 2 = J{I- I a ) dt, 

t r tl 

t t t- +00 

«. = /(./- 4)*, « = -Ji.«-Jt.« + eK- < 

t— — oo 
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We have the identity: 

I + + eu + = (J+) = J r + J jdt 



tr 

t+ u 



I r + eu r + / ( j - j a ) dt + / j a dt 



I r + eu r + e 1 + J (i a + eu a ) dt 



t+ 



I r + eu r + e 1 + / I a dt + eu" - eu a r 



t r t + 

Ii + J I dt + J i a dt + eu r + eu\ - eu a r + d 

t[ t r 

Ii +em + J (I - I a ) dt + J I a dt + J i a dt + eu r - em + eu\ - eu a r + a 

tl tl t r 

h t T t+ 

J_ + £tt_ + / jdt + / I a dt + / I a dt + eu r — eu\ + eu\ — eu a + e 1 + e 2 



t— tl t r 

tl t r t + 



Z_ + eu_ + / J a dt + / I a dt + / I a dt + en r - eu; + en" - eu" + e x + e 2 + e 3 



t+ 

J_ + + J i a dt + EU\ — + £W r — + — eu a r + ei + e 2 + e 3 
t_ 

+00 

I_ + eu_ + / J a dt + eu r — eui + eu? — eu a r + e± + e 2 + e3 + e4. 



-00 



333 
We should prove that eu r — eui+euf— eu^ = 0{£2), ei+e 3 = 0{£2), and e.,- = j = 2,4. 
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4.2.2 Estimate of eu r — eui + euf — eu^ 

f) **1 dS a 

Here we use the sum u = U\ + u 2 + u 3 + w 4 , u a = u± + u\ + u\ + W4, = — -g- 4 -, = — g-S 
k = 1,2,3,4. For convenience, we consider S\(J,<f,r) as the main part of expression of 
S(J,(p,r), then consider the others. 



1 *-„.a, 



eui r — eun + eu ll — eu lr 

dS lr dS u dSf, dS? 
-e — — + e — — - e — — + e — — 

dp dp dp a dp a 

u)(r r ) uj(ti) Vl{ti) Q(r r ) 

^(r^Q^ J ; , pi, ti) - ^(r^Q 1 ^, (p ai ,n) _ ^l(r r )a l (J r , <p r , r r ) - u(r r )a l (h, p ar , r* 



u(ti)Q,(ti) w(r r )fi(r r ) 
_ [u(ti) - \u"{ti - t*) 2 + 0{n - T^fja^J^ip^Tj) - u(T{)a l {h,p an T,) 

u{ti)VL{ti) 

[u(r r ) - \u)'l(r r - n) 2 + 0(r r - r*) 3 ] a l ( J r , tp r , r r ) - u^a 1 ^, <p ar , r* 

u;(T r )f2(r r .) 



— £- 



w(Ti)f2(r;) 

T r ) [a 1 { J r , ip r , T r ) - Ot l (h, p ar , T*)] - ( 7V ~ T *) a; 1 ^, <£ r , T r ) + 0(t t - T*) 



a l {J r ,p r ,T r ) - ^(I^Pa^T*) Uj'^T r - T*) 2 a\J r , p r , T r ) . 3 

E fi(r r ) 2a;(r r )0(r r ) l£2j 

(a^J/,^,^) - a 1 ^,^,,^)) + (a 1 (J n <£> n T r ) - a x (h, p> ar , t*)) 



£wf(7) - r*) 2 f a l (J u p U Ti) a l (J r ,p r ,T r ) 



0(d 



2 0(n) V w(rj) w(r r ) 
0(VS)E 1+ + 0(sl) f + ^(^^r ; ) + Q(Vi) \ + s 

0(Vi)B 1+ + 0(*l) ( a' W , w , n) ■ + + 0( £ f ) 

V cu(ri)cu(r r ) / 



0(v^)^i+ + 0(e 



Here E i+ = J h p h n) -a 1 (I*, p av r*)] + J r , <^ r , r r ) -a 1 ^*, v? ar , r*)] , and we are going 
to estimate E 1+ using Lemmas 4.2 and 4.3: 
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d^jl^ipa^n) da l (h,Lp an T*) d^jl^tpa^n) 
Ei+ = ^ {J t -Q + — {ip l -i Pai ) + — {n-n) 

da^I^tp^r*) f da 1 (I*, ^.,7-*) da l (h,ip ar ,T*) 

H qJ [Jr - h) H q- [ifr ~ <Pa r ) H 7^ (t p - T*J 

+ 0(J, - J*) 2 + 0(ipi - ip a f + 0(n - r*) 2 + 0(J r - Q 2 + - ip ar f + 0(r r - 

= 7^ [{J l ~ h) + (Jr - ^*)J H - ¥>a,J + (<^r - W 



<9r 

0(e) 



[(•n - r*) + (r r - r*)] + 0(^ aj - p ar f + 0(e) 



Therefore, eu\ T — euu + eu\ { — eu\ T = 0(e*). 

Then we consider the term S 2 (J, <p, r) in expression of S(J, cp, r). Here 

dS 2 _ u'(T)a 2 (J,<p,T) | /3 2 (J,^r) 
<9<^ u 3 (r) Uj2 ( T ) 

dS% _ Q'(r)a 2 (h, Va ,n) 
dp a fi 3 (r) 

eu 2r - Eu 2 i + £u 2l - eu a 2r 

dip dip dip a d(p a 

^ 2 Uj'(T r )a 2 (J r ,ip r ,T r ) _ ^ 2 f3 2 (J r ,ip r ,T r ) ^ 2 u'(Ti)a 2 (J h ip U Ti) ^ 1 (J h ipi, Ti) 

uj 3 (r r ) uJ 2 (r r ) uj 3 (n) u 2 (n) 

^ 2 n'(Ti)a 2 (I*,(p ai ,T*) | ^ 2 tt'(T r )a 2 (I^ip ar ,T*) 



From Lemmas 4.2 and 4.4, 



— £ 



2 P 2 (J r ,ip r ,T r ) 2 f3 2 (J h (pi,Ti) 



+ £' 



CU 2 {T r ) CU 2 (ti) 
2 /3 2 (J r ,ip r ,T r ) 2 f3 2 (J h ip h Ti) 



n 2 (r r 



£ 



n 2 (n) 



[f3 2 (J h tp b n) - f3 2 (J r , <p r , r r )] + 0(e*\ 



'dp 2 dp 2 dp 2 

0(el) 
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Also we know that 



fi 3 (r) = [ w (r)-^:(r-r,) 2 + 0(r-r,) 3 ] 3 

= ^(r) - ?«V(t)(t - + 0(r - r*) 5 . 



So 



2 n 3 (r r )uj'(T r )a 2 (J r ,ip r ,T r ) - uj 3 (T r )tt'(T r )a 2 (I*,Lp ar ,n) 3 
£ io 3 (r r )n 3 (r r ) +V{E2) 



s 2 



[u 3 ^) - §o;V(r t )(r t - r.) 2 + Q(r t - T.)V(r t )a 2 (.7 t , Vl , n ) 

(n) [cu'jri) - u"(n - n) + Q(n - n) 2 ]a?(I*, tp ai ,n) 
co 3 ^ 3 ^) 

2 [u 3 {r r ) - |w> 2 (r r )(r r - r*) 2 + 0(r r - r*) 5 ]u/(r r )a 2 ( J r , <p r , r r ) 



9 w 3 

- £^ 



a; 3 (r r )fi 3 (r r 



2 w 3 (r r )[a/(r r ) -u"(r r - n) + 0(r r - r*) 2 ]^/*, y? ar , r*) 3. 
+ £ 3 , x n3 , x J + 0(e*) 

u 3 (n)u'(n)[a 2 (Ji,ipi,Ti) - a 2 (h,if ai ,n)] 



s 2 



a; 3 (r z )n 3 (r z ) 



e 



u; 3 (t*)0 3 (t,) w 3 (rOn 3 (rO 

2 w 3 (r r )a; / (Tr) [a 2 (</r, <fr, T r) ~ a 2 {h, fa r , T*)] 
w 3 (r r )fi 3 (r r ) 



2 



2 |u;V(r r )a; 2 (r r )(r r - nfa 2 { J r , yy, r r ) 2 ^a; 3 (r r )a 2 (/„ yy, T«)(r r - 7%) 3 

u;3(r r )n 3 (r r ) 5 w 3 (r r )ft 3 (r r ) l£ 

a;' co 1 
e 2 -^-j[a 2 (J h <p h T l ) - a 2 (h,ip an r*)] - £ 2 ^ 3 * . [a 2 (J r ,cp r ,T r ) - a 2 (h,ip ar ,r*)] 

3 2 w''u)'(Ti) . .9 o/ t \ 3 2 w"u'{t t ) . 2 o . 



(ri-T*)a (I*,ip ai ,T*) - £ nv * , (r r - r»)a (I*, w, r*) + 0(^2 ) 



n 3 (r,) v y v ' r " ' ft 3 (r, 
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£ 2 < 



(a?(Ji,ipi,Ti) - a 2 (h,<f ai ,n)) + (a 2 (J r ,ip r ,T r ) - a 2 (h,ip ar ,n)) 



2 ^* ( 7~l ) r 2 / t \ 2 / r \ 1 2 ^* ( 7~* )r2/r \ 2/r 

^3( T; ) ^ W'^' T «) ( J *^op T *)J _£: ^ 3( y j L a { J r^r,T r ) - a (h,ip a , 

3 2 W * W *( r i - r *) 2 r 2/ t \ 2/ t M 

+ g2 ^%7^ [« 2 ( J - ¥>«,. r.) - a 2 (/„ ^, r.)] + 0(e§) 
= 0{V~e)E 2+ + 0(e) [a 2 (J u <p u n ) - a 2 (h, cp ap n)] + 0(e) [a 2 (J r , tf r , T r ) - a (I*, (f ar , T*)l 

+ 0(£)[a 2 (J i ,^,r i ) -a 2 (J r ,^ r ,r r )] + 0(e) [a 2 (i*, ^ ap r*) - a 2 (i*, <^ ar , r*)] +0(^5) 
= 0(^~e)E 2+ + 0{el), 

where £ 2 + = [a 2 (J h ip h n) - a 2 (/*, <p av r»)] + [a 2 ( J r , y2 r , r r ) - a 2 (/*, <p ar , r*)] . 



Similarly to Ei + = 0(e), we can obtain E 2+ = 0(e) from Lemma 4.3. Therefore, it is true 
that eu 2r — £u 2l + eu 2l — eu 2r = 0(e*). 



Similarly, we can derive that 

eu ZAr - eu 3Al + eu% M - eu% Ar = O(ei). 

Therefore, 

3 

eu r — eu\ + eu^ — eu a r = 0(e?). 

4.2.3 Estimate of e 2 

We simply use Lemma 4.8 and Lemma 4.6 in order to get estimate: 

e 2 = J(I-I a )dt = -ej { j dt = 0(e 

ti k 

4.2.4 Estimate of e 1 + e 3 

For combined term e 1 + e 3 , we consider j = j x — j 2 + 0( ( r _ £ J ^ 10 ), where 

^ = ^ 7(^(r))V(J,^,r) + c5 l (J,p,T) 
u 8 (t) u 7 (t) 

j ( 5 7{u'(T)) 4 a\J,<p,T) 5 7 (J,y,r) \ A fc 7 fc (J,y,r) 



19 



and also 

j a = £ 



B 7(fi'(r)) a A (I*,<p a ,T*) 5 7(J*,^ a ,r* 



ft 8 (r) fi 7 (r) ' 

For Ji — J a , we apply Lemma 4.1, then Lemma 4.9: 

Ji-Ja = e [a (J, y, r) - a (J., y tt , r.)] + s ^— ^ - e ^ + O { 

T 

— dri 



7* 



+ 7 « + 7 ^ / dr l 



r - r*) 5 (r - r*) 8 J dl (r - r*) 7 



+ o [j ^)+o [j7 ^ )+ o [j —^ )+ o [j — 7]+ o 

e 6 5 5 (h,ip a ,T*) e 7 5 6 (h,ip a ,T*) e s 5 7 (h,ip a ,T*) e 9 5s(h,(p a ,n) 
(r-r*) 9 + (r-r*) 11 + (r - r*) 13 + (r-r*) 15 

Here 7 is a smooth function. Also 

. / 5 7Hr))V(J,y,r) . 5 7(^,t) \ A r) 



fc=l 



fc = l 

4 



I r T - r J 2fc + 7 ^ (r- rj 2fc + 6 7 Ifr - rj 2fc + 8 7 V fr - rj 2fc + 4 



y- fc+5 7 fe Cr*> ¥>q, r,) y> / £ k+i \ 
ti (T-^) 2fe+7 V(r-r,) 2 ^; 



Using Lemma 4.5 for estimate of integrals and Lemma 4.7 for cancellation in symetric inter 
vals, we obtain: 



t+ t t t+ ti 



ei + e 3 = y ( Ji - j G ) dt + y ( jj - j n ) dt - y j 2 dt- J J 2 dt + 0(el) = 0(e 



tr 
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4.2.5 Estimate of e 4 

For term e 4 , we apply an integration by parts: 

- /~/ n d* = e *j ^(1 <p a ,r.) df 



dtp 

— oo 

t- „ 

f dH^I^ip^T*) . 1 

E J * v '-W) At 

— oo 
t- 



eH 1 {I*,(p a _,r*) t ^ 2 /" Hi(I*,y a ,n) ^ 



— oo 

eHi(I*,<p a _,n) 



0(r_) 



and similarly 



+oo 



+ 



Therefore, the estimate e 4 = 0(e 2 ) is obtained. 

By joining estimate of terms ej together, and taking into account that u± = u\± + 0(e), as 
well as the identity 

+oo +oo +oo , 



dip J dip 

— oo — oo — oo 

where # = — =— , we have finished the proof of the first formula of Theorem 2.1. 
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4.2.6 Principal identity for formula (2.6) 

Let 

u u t- 

ei = J{I- /«) dt, e 2 = J (J- j a ) dt, e 3 = - J I a dt. 



We have 

u 

h = Ii + J idt 

tl 



/i + J{i-i a )dt + J i a dt 

ti ti 
i, 

Ji -eu t + J I a dt + ei 

J_ + / J dt - euj + / J a dt + e x 



t- h 

J_ +EU- + J(j~ -Ja) dt + J J a dt- eui + J I a &t + ei 

t- t- ti 

I- + J (i a + eii a ) dt + eu^ -eu t + J I a dt + ex + e 2 

t- ti 

I- + J I a dt + euf — eu°i + eu_ — eui + §i + e 2 



7_ + y I a dt — J I a dt + en" — en" + — eu t + ei + e 2 

— oo — oo 

i* 

I -e [ dHl ^*' ^n) dt + eu? - eu a _ + - em + ex + e 2 + h 
J dip 

— oo 

+oo , 

T 1 r f dHx{h,y* + ^6 2 ,T*) jn a „ _ _ 

I- - -ye / o d6> + eu, - en_ + eit_ - eu; + e x + e 2 + e 3 . 

2 y ay? 

— oo 
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Making use of Lemmas 4.1, 4.2, 4.4, 4.5, 4.8, 4.9, one can show that e.j = 0(e), j = 1,2,3, 
eu? — eui = 0(e), eu- = 0(e), eu a _ = 0(e). This leads to the first formula of Theorem 2.2 
with the sign "— ". The proof of the formula with the sign "+" is completely analogous. 



4.3 Proof of formulas for angle variable 



4.3.1 Principal lemmas 



Lemma 4.10. For fcGN and smooth function f(J, tp, r), 



e k f(J,V,r) 



(b) 



U 2k-2( T ) 



e k f(J,^,r) 



uj 2k -Hr) 



0(e 



e k f(h,(f a ,T* 



T 



+ 0(e*). 



Proof. For r E {r h r r }, Lemma 4.9 can be simplified as f(J, ip, r) = /(/*, <p a , t*) + 0(r - 
Thus with Lemmas 4.1 and 4.4, 



(a) 



s k f(J,^,r) 



uj 2k - 2 (r) 



e k f(J,V,r) 



Q 2k ~ 2 (r) 
e k f{J,V,T) 



e k O(l) 



n 2k - 2 (r) 



Q 2k ^(r) 
+ 0(el) 



£ k f(I*,<fa,n) 



Q 2k ~ 2 \ 



T 



tr 



e k O(l) 
Q 2fc - 3 (r) 



0(e 



0(e 



Applying the result of (a) and Lemma 4.9, we obtain 



(b) 



uj 2k -Hr) 



e k fU<P,r) 



e k f(J,^,r) 



e k f(J,^r) 
n 2k ' 2 (r) 



+ 



e k O(l) 
Q 2k - 3 (T) 



n 2k -h 



T 



0(e 



£ k f(l*,<p a ,n 
e k f(h,ip a ,T* 



e k f(h,(p a ,n 
n 2k ~ 2 (r) 



0(e 



+ 0{e*). 
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Here / and / is smooth functions. □ 

Lemma 4.11. Let f = f(J,<p,r) be a twice continuously differentiable function. Then for 
k E N, k > 5, 

£ k f(J, <P, t) ^ , f e k f(J, <p, r) ^ _ ^3 



t- tr 

} e k f(J,<p,T) J e k f(J,<p,r) J e k f(h,<p a ,n) f e k f(h,<p a ,n) j 

{b) J Z^=Hr) dt + J Z^=Hf) dt ~J n^(r) dt + J fi»-2( T ) dt + 0{£ ' 

Proof. For t G [£_,£/] U [£ r ,£ + ], we can simplify Lemma 4.9 as follows: 

f(J, V, r) = /(/*, <p a , r,) + 0(r - n) + O ^ 'J*^ + O 
Thus with Lemmas 4.1, 4.5 and 4.7, 

fa) /' efc/ i J ^\ r) dt + / ^44^ dt 



w 2fc - 3 (r) J w 2fc " 3 (rj 

V/(j,p,t) / ^,y,r) ? e fc Q(i) d / ^Q(i) d , 

tl t+ 

e k f(J,ip,r) , /" e k f(J,ip,r) , s, 

V/<w.) ?^ow d(+ /£^ d(+ /£!?m d( 



t_ t_ t_ t_ 

+ liEM^ At+ '(^L At+ 'f^M + Jgrnv + otf- 

J Q 2k ~ 3 (r) J fi 2fc ~ 4 (r) y fi 2fc ^ 6 (r) y fi 2fe - 9 (r) v ' 

0(el). 
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Applying the result of (a) and Lemma 4.9, we obtain 

(b) t [ £kf ^ r) at+ ? £fc/(J ^ r) dt 

W J u^(t) m + J u^ 2 (t) m 



t- t- t_ 

7 e k f(J,(p,r) , /" e k f(J,ip,T) , */ e fc O(l) , 
7 fi 2fc - 2 (r) J fi 2fc ^ 3 (r) J fi 2fc ~ 4 (r) 

£y £7" 

e k f(J,ip,T) , [ e k f(J,ip,T) , ~, 3, 

i#^) i d'+/ip^r di+0(£5) 



t_ t_ t_ 



2fc - 2 (r) 7 Q 2fc - 3 (r) 7 fi 2fc - 5 (r 

ij" 1 

+ e*- 2 0(l) 



Here /, / and / are smooth functions. 
Lemma 4.12. Fort G [t-,tj], 



i_ 

Fort G [t r ,t+]j similarly we have 



■ 



t 
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Proof. Making use of of formula (3.18) and Lemma 4.9, for t G [£_,£;], we get 

t- t- t- 

t 

_ 5 7o;^ a 4 {I*,<p a ,T*) ( e 4 



g — n 7 ; d ^i + 



^ 8 (ri) 1 V(^-^) 6 

Similarly for i G [t r ,i+], we obtain the second formula. 

4.3.2 Principal identity for formula (2.5) 

The relation between tp + and is: 

t+ 

<f + + £V + = (lf) + ) = 1p r + Jlpdt 

t+ 

= (f r + ev r + ipdt 



tf tf t-^- 

<Pi + J w(t) dt + J e^El^pil dt + ev r + J ip dt 

t% t% t r 

4i l -evi + ev r + J oj{t) dt + j £ ^liL^flll dt + J ipdt 

(f^+ev^+ J ipdt- sin + ev r + j uj{t) dt + j £ ^EAL^lll dt + j ip dt 



t- t t h 

t+ t r 



oj(t) dt + / e- 



9/ 



dt 



t; ti 



t- t_ 



k=5 
9 
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Let 



<9#i(/,yy 
dl 



k 
k 



k *+ 
dt, ej +1 = J ^-/K-i dt + y A^'ft,- dt, j = 1, 2, 3, 4, 



tl 9 * + 



fc=5 



The identity becomes 



ip + + ew+ = f- + £i>_ + - / ^( r ) dr - ew/ + ew r + ef + e% + e% + e\ + e% + e% + e£. (4.27) 



We will discuss the estimate term by term. 



4.3.3 Estimate of —evi + ev r 



Here v = V\ + v 2 + vs + v±, = k = 1,2,3,4. Making use of equations (3.17) and 
Lemma 4.10, we get 

—evi + ev r = 



dSi dS r 



d,J dJ 



e dH 1 



ui dJ 



J dH l 



e 3 ap 3 



+ 



u 4 dJ 



uj 3 dJ 



,3u/ 2 dH l 



dJ 



— e 



15u 13 dHx 



cu 7 dJ 



e 2 d(3 2 



+ 



u 2 dJ 



+ 



dJ 



e dHi(I*,<p a ,T» 



n di 



fl' dHi(I*,<p a ,T* 



O 3 dl 



,3ft' 2 dH^p^n) 



n 5 



dl 



— e 



15ft' 3 dH^h^r*) 



ft 7 



dl 



+ 0{£*). 
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4.3.4 Estimate of e 



Using Lemmas 4.6 and 4.8, we have 



^4 
ol 



tt 



dl J dl 

ti ti 

J dl J dl J dl 

ti t* ti 

£ J di dt + £ J di dt + £ J 97 dt 

t, *» ti 

)d 2 #i(J_,r) /r T ,, */ d 2 H l U+,T) fT t \ ~\ 3, 
+ e y ^-^(J- J_)d* + e y ^±^(/_ J + )dt + 0( £ 5) 



dl J dl J dl 



+ e 
ti 

Lemma 4.13. 



/ -J-)*t + eJ **^(/ - J + ) dt + 0( £ i). (4.28) 



£ J(I-Jj)dt + e J (I- J+) dt = O(et). 



t. 
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Proof. Making use of Lemmas 4.8, 4.9 and 4.12, we get 



e / (J - J_)dt + e / (/- J+)dt 



e J{I-I l + I l -J l + J l -J.)6t + e J(I-I r + I r - J r + J r ~ J+) dt 



e / -e 



t r / t 



1 dHijl^n) dSi 
dti + e— - 

dp dip 

t 



Jdt dt 



+ e 



dJh{l^n) dS r 
e i dti + s- 



dip 



dip 



Jdt dt 



dHi(h,ip a ,r,) a l {J h ip h Ti) 2 u'(T l )a 2 (J h ip h Ti) 

— utl + £ ; — ; h £ 



k L h 



dip 



(u'(ti)) a 3 (J h ip h Ti) 4 (u'(ti)) a 4 (J h ipi,Ti) 



+ £ 



Jdt 



dt 



— uti + £ — ; r £ 



+ £ 



dip <jj(t t ) 
t 

3 (u'{T r )) 2 a 3 {J r ,<p r ,T r ) 4 (u/(r r )) 3 a 4 (J r ,y2 r ,T, 



o) 5 (r r ) 



LU'[T r 



cu 3 (r r ) 



+ J , 



dt + 0(e 



^(I^ipa^r*) ^ 2 ^a 2 (h,ip ai ,T^) a 3 (h,<p ai ,n) ^Vt'^ a 4 (h, ip ai , r*) 



ti 



_ 5 7^ 4 a 4 (J*,v2a,T-* 



fi 8 (T!) 



dti + £ 



n 3 (r r 



+ £ 



3 tt'^a 3 (I*,ip ar ,T*) , ^4^ 3 a 4 (/ !l ,,v?a r ,r* 



n 5 (r r 



+ £ H 



fl 7 (r r 



Ju'^a^jl^ip^T* 



dti 



dt + 0(e 



□ 
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We also have (cf. Section 4.2.5): 



tr 



+00 



01 
01 



dt _ £ d Hl (I y tt) n) dt _ £ f dH x {h<p» t., dt 



01 



01 



t~ „ +00 „ 



01 



01 



+00 ~ ~ *+ ~ 



01 



01 



01 



e dHi(h,<p a _,T*) e dH^I*, <p a+ , t*) 



+00 



01 



n(r+) 



9/ 



'+ 



+ 



3/ 

£ &Hi(I»,<p a ,T*) 



n(r) a/ 



d#-£ 



0(. 2 ). 



a/ 



dt-e 



9/ 



dt 



Using this result, formula (4.28) and Lemma 4.13, we obtain 



>''^ T) dt +£ fo-<V>dt 



+00 



9/ 



9/ 



97 



d# 



01 



01 



n 01 



1+ 



0{e 
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4.3.5 Estimate of 

We know that TZi = (Hif = H\. Thus with Lemma 4.12, we have 

k t+ 

J dJ J dJ 

t_ t r 



dt + £ J — ^p^( J - J -) dt + e J o(j-j.ydt 

t- t- t- 

t+ [ dHAJ + ,T) , fd 2 H 1 (J + ,r) fT TN1 V , T 

+ £ / — kr dt + £ / — g/2 ; (J-J+)dt + g / o(j- j+f&t 



dl J BP J n 8 (n) 

t- t- t- 



4.3.6 Estimate of e£ 



From Lemma 4.9, we get the error between J and 



<9y9 r — r* (r — r*) 3 (r — r*) 



5 



+ e 4 J ) * ,ra I* 7 + 0(e) + 0(r - r*) 2 + O ^ ^- + O 



;t — r*) 7 \ e j \ e 2 



T T* 
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As 



dH x dSA* _/ dH l a l \ v _ I dH 1 ( -E x 



dl d<p / = \ dl u / \ dl 
I/ 1 d fA*- 1 9 



we have 



s 2 



!aj <9/ 2 a +£ J V 2cJ dP a 



2uo 

t- t 



+£2 ?/n^y„r) ttt , /f_M^m.r) +)dt+ ^ 



2^7 a/ 2 J \ 2u) dP 



2wy ar 2 J \ 2loJ dp J n 8 (n) 

t- t- t- 

+ 0(e 2 ) 



2^7 &r 2 ~" 1 ~ y v 2^7 &r 2 
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p.v. e 



t- 



aJ <5->o+ 



1 \ &(Hiy{J_,n) 

2uJ <9/ 3 



t*-<5 



1 \ d\Hiy 



2cu J dP 



dt + e z 



1 \ d 2 (HlY 



2u / dP 



dt 



U+5 



+oo 



i a dt + 0(e*) dt + 0(e 



tr 



t-L . — . +00 ~t-\- 

p.v. e I — dt + e 



dl 



op 



dip 



dt / — dt + 0(e 

Id 



+oo 



p.V. £ 



dl 



dt 



4w' 



dP 



dip 



dt + 0{e*). 



4.3.7 Estimate of el 



dH l dS 2 Id 2 ^ fdSi 



2rr , a n X 2\ V 



So 



dl dip 2 dP \dip 

dHi-u/Hi 1Q(1)\ 
~dl Z^ + 2~^J r J 

ur \ dl cu 2 



0(e 
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4.3.8 Estimate of e£ 



, dH x dS 3 , d 2 H l dSx dS 2 , 1 d 3 H l fdS^ 3 

/V4 



dl dtp dP dip dip 6 dP \ dip 

dH 1 -3aAgi 0(1) 1 0(1) 
dl co 5 to 4 6 to 3 

3oo' 2 /^dH i y 0(1) 



cj 5 \ 9/ / u 



So 



, 3oo' 2 d /*~dH x y , T , , 4 /-3a;' 2 d / £ T dH l \ v , T . , 3, 
-e 4 (H,— i) (J,r)dt-e 4 / — — (#!— -i) (J, r) dt + 0(e*) 



w 5 <9J \ <9I / v ' ' J u 5 dJ \ dl 

,}m ,2 d/%dH 1 Y fT NJ 4 t f3Q' 2 d/^ T dH l \ ip /T NJ ^, 

0(e§). 



4.3.9 Estimate of e£ 

There exists a function [/(/*, </? , r*), such that = — — 1 ^ *' *^ . Here 



C/ = ^Elil^pllll with = 0. Thus 



dH^h^r*) fdU fdU 1 /"I 

£/ 57 dt = £ J d^ dt = £ J ^-n dt = £ J n dU 

eU{P,ip a ,n) t ^ 2 fjj^ dt 



fi(r) y fi 2 

There exists a function V(J*, </? a , r*), such that ^— = U[I*,<p a ,T*). Here 
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V = U(h, <p a , r.) = gglfejg^j with {V y« = o. Thus 



TT Q' 2 fdVQ' 2 fdVn' , 2 /" fi/ , T , 



There exists a function Wm/*, r*), such that — — = V(i*, u> a ,T*). Here 
W = t>(/*,^,r*) = d Jh^LpllA with = o. Thus 



."/V-jjj-d* = *J——« = *J——« = Sj—iW 



dT 

There exists a function T(l*, y? a , r*), such that = <p a , t*). Here 

T = W(h, <p ai r.) = ^lil^pl^l with (T) ^ = o. Thus 



= ^'••™)-w?; L + <?J T -v- d( - 

Thus we have 

r eE l{ i.v.,T.) a = eu + £2v v +e3w zv> +elT ^ + eS r T m^ dt 

J 81 Q fi 3 fi 5 fi 7 7 ™ 

5.9^1 ~VL' dH 1 o3fi ,2 9#i .ISO' 3 5^ - /" 105^' 4 
= -+£ -+£ -+£ - + £ / 

(is/ di n 5 di n 7 di J n 8 di 
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Therefore, applying equation (3.16) and Lemma 4.11, we obtain 



.dS 4 , / d K dS, 



s 5 ^dt 



dJ dr 



t- 



k 5 105cy 4 djw^r) , } s 5 /3(J,ip,r) 
£ n QZ - 



dJ 



7 



dt 



dJ 



e 5 P(J,<p,r) 



dt 



£ / — — ^-r^ i dt + 5 / — — —rz -dt + 0(e 



n 8 



dl 



f dHi{h,y a ,n) 
dl 



e dHi(h,y a ,T* 



n 8 di 

t, 



dl 



s 



t- 



01 



— e 



+ e 



dl 
dl 



t- 



dl 



dt 



0(r) 



dl 



01 



,3fi /2 (r) dH x {L,<p a ,n 



dl 



— e 



15fT 3 (r) dH 1 (h, l p a ,r, 



dl 



dl 



+ 0(52) 



9/ 



n di 



fi' dHi(I*,<p a ,n 



n 3 di 



_ 3 3fi' 2 dHi(I*,Lp a ,n) 



n 5 



di 



15Q' 3 dH x {h^ a ,n) 



n 7 



01 



0{e 



where /3 and /3 are smooth functions. 
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4.3.10 Estimate of e£ 



Applying Lemma 4.11, we easily get the estimate 

Sr ^ k U^r) f ^ k 7 k (J,<p,r) i 



Combining obtained estimates and the identity (4.27), and taking into account that v± 
v\± + 0(e), we get the second formula in Theorem 2.1. 



4.3.11 Principal identity for formula (2.7) 

With results of four steps of perturbation theory, we have 



v. = <p l+ MT)dt+ e ; dHl ^ T) dt 



dl 

k k 
k t» U 



I ip dt - evi + / u(t) dt + e I dt 



t- h t t 

U dH^I^r) 



+ £u_ — ev i + y w(r) dt + e J 
t- h 



9/ 



df 



+ / — I + e'^ 2 + £^ 3 + £ 4 ^ 4 + ) dt + 0(e) 



t- 



From this relation, similarly to the above proof, we obtain the second formula of Theorem 
2.2. 
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5 Numerical verification 



5.1 Example 

We will verify previous formulas numerically using the following example suggested in [3]: 

(5.29) 



< dJ Mr ^ 

— = -eA(I,t) cosp 



. . dA(I,r) 



where 

w(r) = e T_1 - 1, A(I,t 



Ve 7 "" 1 + 1 



We choose slow time interval r G [0,2], take 1(0) = 1 and consider different values <p(Q) G 
[0,2tt]. We have 

H^I, p, t) = A{I, t) sin tp, H x = (H x y = 0, H X = H X -E X = H x , 
Hi A(I,t) sin if> du\ dA sin if dv\ 



Ul 



cu e^" 1 - 1 ' dl " dl e^" 1 - 1 " dp 

I 



and 



OA cos if dA(I,r) _ 2\/4 — i" 

Vl ~ We^ 1 - 1' 97 ~ Ve"" 1 + 1 

i_ = 1, r_ = 0, r* = 1, r+ = 2, cu* = 1. 



5.2 Theoretical value of I + 

According to (2.8), we should find 



Tto** = i_+e Ul (I_,<p_,T-)-eui(I-,<p + ,T + )-V£ / tUHZl &e (5>30) 



dHi(h, p* + ^-6 2 ,t* 



df 



We perform as follows: 
1°. Value of ip+. 



T+ 2 

If... If,,, e -e" 1 - 2 



<p + = <p- + - J lo(t) dr = <p- + - J oj(t) dr = tp_ + 



o 



3s 



2°. Value of <p*. 



T„ 1 
If I f 

- / oj(t) dr = <f_ H — / w(r) dr 



y?_ . 



3°. Value of </?*. 

1 ? , 1 r + f > dH 1 {I.,^ + ^B\r,) AQ £ ln £ a 2 (^y(J-,< 
y?* = + - / w(r) dr + - ve / — — d# + 



dl 



4uL 



dp 



Here 



+oo — , 

2^ J dl d ° 

— oo 

+oo ^4 _ j L — 

— oo 



V2 



+oo 



cos sin y + sin 0* cos y 



d0 



and 



2\/6 
5^ 



4\/6 

ar 2 



+oo +oo 

cos 0* y sin y d^ + sin ^ cos y d# 

o o 

cosy5* + sin^*). 



d 2 /I 2 (A-I) 



dl 2 \ e^" 1 + 1 

i d 2 / 2 (4 - /) 



sin 2 y?) (i_,T*) 



2 &T 2 e T - x + 1 



1 8-6J_ 1 



2e r *- 1 + l 2' 



Therefore, 



<p* = tp* -\ 7=r- (cos + sin y?*) + 



4^ 
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4°. Value of I*. 



+00 / 

* T i r f diui^+^e 2 ^) 
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Here 



+00 



1 f dH 1 (J zl <p 1 + ^9 2 ,n 
dip 



2 V£ 



d0 



-00 
+00 



/ 7 V4 1 1 cos + ±6 2 ) &6 



+00 



V3 



2a/2 



e 2 ■ ~ ■ e 2 
cos </?* cos -2 — sin sin y 

cos 0* ■ - a/2 sin • 
cos (,5* - sin^*). 



d9 



Therefore, 



We have also 



i* = 1 ^=r- (cosy?* — sm^*, 

2 v 2 



, T v v^£sin^_ , r „ . v / 3esin(^_ + ^^) 

(e _1 — l)ve _1 + 1 (e — l)ve + 1 

and 

+ OO v , +OO v , — 

V£ / o = v £ / — , : cos + -^-0 I dg 

J dip J Ve T * _1 + 1 



-00 




Thus we have formulas for all required values for calculation of /™ cor (5.30). 
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5.3 Theoretical value of ip + 



vT m = ¥>- + eviV-,<P-, r-) ~ r + ) + -f w(r) dr + y/i [ dH ^<P* + e 2 ^) ^ 



From (2.9), taking into account that Hi — 0, u'^ — 1, we get 

T+ +00 

dHi(h, 

dl 

T— — OO 

+ p - v - £ i aj dr — i op J d~ v de - (5 ' 31) 

r_ — oo 

We perform as follows: 
1°. 

1 /" e — e — 2 

V?- + - / u;(r) dr = 0+ = y?_ H . 



2°. We know that 

cL4 cos</9 / , / \ cosy? 



Thus, 



" 81 e- 1 - 1 V 2 v / 4^T; (e^ 1 - l)y/&-i + i 

£v(I T ) = (^-^k) £cos ^- = i75 £CQS( ^- 
evi{ ,<p , r_j ^ _ 1)v / e _ 1 + 1 ^ _ 1)v / e -i + 1 ' 

v (V3-^=)ecos0 + ^j£cos0 + 



3°. 



(e - l)v^TT ' (e - l)Ve~TT 



+oo ~ ~ +00 . I A J _ I* 

e / m d6 = ^ej v ^-rr TT sm ^* + T ) dg 



— oo 



8 — 3 J* „ /- /7F, v . „ N a/7tF(8 — 3/*) . 

• ; • 2V2 ■ \ —(cos if* + siii(/3j = (cosy?* + sini/3, 

L V 8 



2V2 V4 - J, V 8 V8 V 4 - 1. 

Values /*,</?* were calculated in the previous subsection. 

4°. We have TZ 2 = - — ^-(HfY = - — 4 ( ~ sin V V- Therefore 

<9ft 2 1 d 2 / ~ 2 , ip _ 1 8-6/1 1 4-3/ 



9/ 2udP x 17 2w e^ + l 2 2u e T ' 1 + 1 
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Hence, 



I Qj dr 



l \ 4 " 3/ - A 

p.v. e I I : dr 

2uj ) e— 1 + 1 



P ' V - " 2 / e— 



- dr = p.v. 

1 y 2 



-:2t 



q2t _ e 2 



1 dr 



■— lim 

4 <5^0H 

■— lim In 

4 <5^0+ 



|e 2r 


-e 2 | 


1-5 

+ In |e 2r - 




-e 2 | 


2 

l+<5 






-e 2 


— — lim In 

4 s^o+ 


e 4 - 


e 2 












+ £ 


,2+2<5 


-e 2 


1 - 


e 2 



£ 

2' 



5°. 



Thus 



d 3 ^}* _ d 3 /I 2 (4-J) 



<9 3 J 2 (4 - J) 1 



dP e 



r-l 



1 2 e 



r-l 



+oo 



dP 



2 ' *J 



<9</? 



d9 



4 



— oo 



3v 37r 3 . 
-e 2 me (cos (,5* — sin 



£2 lne 



£2 lne 



8^ 



Therefore, from (5.31), we obtain the estimated value of <£>+: 



e _ e -i_2 . Aecosp_ -5=£cos</?+ 



^ + hcor = P- + — + 



(e- 1 - l)v / e- 1 + 1 (e-l)VeTT 

v / 7rF(8 — 3/*) . v . v . £ 3a/37t s, v . „ 

H . (cosy?* + sin^J + - H ^£2 lne (cos<z>* - sin y?J. (5.32) 

7 4 - L 2 8V2 
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5.4 Results of numerical simulation 



Our goal here is to check numerically, if the obtained values <y2+ eor , /^ heor indeed approximate 
actual values of /, cp on solutions of system (5.29) with the accuracy 0(e^), as it is predicted 
by Corollary 2.1. To this end we introduce the variable 

T 

f W ( Tl ) A 1 f T-1 -1 \ 

X = V ~ / i~ 1 = ^ _ £ ^ ~ G 



and integrate numerically by 4th order Runge-Kutta algorithm the system 

/ = -eA(I, t) cos (x 



e^ 1 - e~ l - t 



£ 

8A(I,t) ( e^-e" 1 - 
X = £ — gj — sin ( x H : 



(5.33) 



on the slow time interval < r < 2. At r = r + = 2 we have values J^ umer ; ^^ umcr and 
^numer _ ^,numer _|_ ^ e _ e -i _ 2)/e. These calculations are performed for 48 values of </?_ 
equally spaced on [0,2tt], and 11 values of e, {0.02, 0.015, 0.01, 0.007, 0.005, 0.003, 0.002, 
0.0015, 0.001, 0.0007, 0.0005}. We calculate 

E I+ {e)= max | J™ r - J| heor |, E v+ (e) = max r - ^ + hcor | , 

i^_g [0,27r] </3_6[0,2-7r] 

and plot values ln(£ , / + ) and ln(£^ + ) as functions of Ine in Figures 1 and 2. Linear least 
squares fit of the data in Figures 1 and 2 gives slopes a = 1.574 and (3 = 1.38, respectively. 

3 

The ideal results for accuracy 0{e^) would be a = 3/2 and (3 = 3/2. Thus the numerical 
simulation indicates that the accuracy is 0(e^), as expected. 
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